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Summary. Complex systems comprise a large number of interacting elements, whose dy- 
namics is not always a priori known. In these cases - in order to uncover their key features 
- we have to turn to empirical methods, one of which was recently introduced by Menezes 
and Barabasi. It is based on the observation that for the activity fi(t) of the constituents there 
is a power law relationship between the standard deviation and the mean value: (Ji oc {fi) a . 
For stock market trading activity (traded value), good scaling over 5 orders of magnitude with 
the exponent a = 0.72 was observed. The origin of this non-trivial scaling can be traced 
back to a proportionality between the rate of trades (N) and their mean sizes (V). One finds 
(V) oc (iV) for the ~ 1000 largest companies of New York Stock Exchange. Model inde- 
pendent calculations show that these two types of scaling can be mapped onto each other, with 
an agreement between the error bars. Finally, there is a continuous increase in a if we look at 
fluctuations on an increasing time scale up to 20 days. 
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1 Introduction 

Although there is no generally recognized definition of complex systems, one of their 
widely accepted properties is that they comprise a large number of interacting con- 
stituents (or nodes) whose collective behavior forms spatial and/or temporal struc- 
tures. Some of them are labeled "physical" because they are treated in the regular 
framework of physics. Nevertheless, the above scheme itself applies to a much wider 
range of systems, including the world economy consisting of companies that trade 
and compete. Human agents can interact with each other, e.g., by social networks 
or on the trading floor. We have little or no a priori knowledge about the laws gov- 
erning these systems. Thus, very often our approach must be empirical. Recently, an 
increasing number of such systems have become possible to monitor through mul- 
tichannel measurements. These offer the possibility to record and characterize the 
simultaneous time dependent behavior of many of the constituents. On the ground of 
these new datasets, an emerging technique (de Menezes and Barabasi 2004a) seems 
to be able to grasp important features of the internal dynamics in a model indepen- 
dent framework. 



2 Scaling of fluctuations in complex systems 



The method is based on a scaling relation that is observed for a growing range of sys- 
tems: The standard deviation <Ji and time average (/*) of the signal fi(t) capturing 
the time dependent activity of elements i = 1, . . . , N follows the power law 

a% oc (fi) a , (1) 

where we define 

- yj((fi-(fi)) 2 ), (2) 

and (•} denotes time averaging. 

This relationship is not unmotivated from a physicist's point of view. The con- 
stant a - while not a universal exponent in the traditional sense - is indeed the finger- 
print of the microscopic dynamics of the system. Applications range from Internet 
traffic through river networks to econophysics. The latest advances (Menezes and 
Barabasi 2004b, Eisler and Kertesz 2005) have shown several possible scenarios 
leading to various scaling exponents: 

1 . The value a = 1 always prevails in the presence of a dominant external driving 
force. An example is web page visitation statistics. Here the main contribution 
to fluctuations comes from the fluctuating number of users surfing the web: a 
factor that is not intrinsic in the structure of the network. The situation is very 
similar for networks of roads or rivers. 

2. There are systems, where the different mean activity of constituents comes ex- 
clusively from a different mean number of events. Individual events have the 
same mean contribution (impact) to a node's activity, only for more active nodes 
more of these events occur. When the central limit theorem is applicable to the 
events, a = 1/2. This behavior was observed for the logical elements of a com- 
puter chip and the data traffic of Internet routers. 

3. Two mechanisms have been documented so far that can give rise to an interme- 
diate value 1/2 < a < 1: 

a) Because of the competition of external driving and internal fluctuations, it 
is possible that cr's measured for finite systems display a crossover between 
a = 1/2 and a = 1 at a certain node strength (/). Then there exists an 
effective, intermediate value of a, but actual scaling breaks down. 

b) The other possibility is related to a very distinct character of internal dynam- 
ics: when elements with higher activity do not only experience more events, 
but those are also of larger impact. We call this property impact inhomo- 
geneity. Stock market trading belongs to this third group with a w 0.72 for 
short time scales (see also Eisler et al. 2005). 

In a recent model (Eisler and Kertesz 2005), the effect of impact inhomogeneity 
has been studied. Tokens are deposited on a Barabasi-Albert network (Albert and 
Barabasi 2002) and they are allowed to jump from node to node in every time step. 
Activity is generated when they arrive to a site. Every token that steps to a node i 



generates an impact Vi whose mean depends on the node degree kf. (Vi) oc kf[ . This 
gives rise to a scaling relation: 

(V) oc (N { f . (3) 
The result of Eisler and Kertesz (2005) can then then be generalized as 



+ 1 

Simulation results shown in Fig^a) are in perfect agreement with formula (|4}. This 
is an example that the value of a is basically determined by this impact inhomogene- 
ity. If (3 — 0, i.e., the mean impact generated on all nodes is equal regardless of their 
degree, one recovers a = 1/2. When (3 > 0, the events on more frequently visited 
nodes are also larger on average. Correspondingly, a > 1/2. 



3 Application to stock market data 

Let us now turn to the case of the stock market. Data for the period 2000-2002 
was taken from the TAQ Database (New York Stock Exchange 2003). We define the 
activity /$ (t) of stock i as the capital flow in time windows of size At. In window 
t, fi(t) is the sum of Ni(t) trading events. If we denote the value exchanged in the 
n'th trade of time window t by Vi(t; n), then the total traded value of stock i is 

/;(*)= X)^(t;n). ( 5 ) 

n=l 

Then, (V) is the mean value per trade, while (N) is the mean rate of trades. 

As we wish to calculate the mean and the standard deviation of this activity, it is 
essential that these quantities at least exist. Traded volumes and consequently traded 
values fi(t) are often considered to have a power law tail (Prob(f > x) oc x~ x ) 
with an exponent Aj ~ 1.5 — 1.7 (Gopikrishnan et al. 2000). This would imply, that 
the standard deviation is already divergent. Recent measurements, however, indicate 
that both of these quantities exist and that there is no unique \ for a stock (Eisler 
and Kertesz unpublished). 

Then, it is possible to test the scaling relation Q and one finds good scaling over 
more than 5 orders of magnitude in (/) with a w 0.72. This is a value which can be 
- at least partly - explained in terms of impact inhomogeneity. We found 1 that for the 
stocks of the - 1000 largest companies of NYSE, (3 = 0.69 ± 0.09 (see Fig.[TJb)). 
Substituting this into ® we expect a = 0.71 ±0.01, which is very close to the actual 
result seen from Fig.|2ja). Note that although large error bars prevent us from testing 
@ for smaller stocks, we still find that the scaling law Q holds. The exponent is 
unchanged, but this can only be explained by a detailed analysis of fluctuations. 

1 The result is qualitatively similar to those of Zumbach (2004) for the FTSE 100. He shows 
that both (TV) and (V) scale as power laws with company capitalization for large compa- 
nies. Capitalization dependence can be eliminated to recover ^3). 




p log <N> (trade/min) 

Fig. 1. (a) The value of a as a function of (3 for the random walk model introduced by Eisler 
and Kertesz (2005). Circles give simulation results, while the solid line corresponds to 0. 
The inset shows actual scaling plots for various values of j3. (b) Plot of mean value per trade 
(V) versus mean rate of trades (iV) for NYSE. For smaller stocks there is no clear tendency. 
For the top ~ 1000 companies, however, there is scaling with an exponent (3 — 0.69 ± 0.09. 




Fig. 2. (a) The scaling a oc {/) Q for fluctuations of traded value at NYSE, At — 5 sec. 
Dots show raw results for each stock (shifted vertically for better visibility), the fitted slope is 
otraw = 0.733 ± 0.004. Diamonds show average cr's for multiple stocks in a range of log (/) . 
This method corrects for bias that comes from the large number of stocks with low (/), one 
finds abinned = 0.720 ± 0.007. (b) The dependence of a on the time window At for the 
NYSE data. One finds that up to At = 1 min, a ~ 0.72, as expected from independent 
approximations. Then by increasing At, the value of a increases. This is due to the presence 
of strong autocorrelations in the activities f(t) stemming from the clustering of trades. 



The mechanism leading to non-trivial a via the scaling can be considered 
dominant only if the events are not strongly correlated. This condition is satisfied 
for short time windows At, when (N) <C 1. Interestingly, the value of a does not 
change noticably up to At ~ 1 min. There is, however, another effect that is rele- 
vant to the value of a for longer time windows. For the NYSE data, a(At) increases 
continuously with At (see Fig. 0b)). Previously (Eisler et al. 2005) this was at- 
tributed to the growing influence of external news: a kind of "driving". With longer 



time given for information to spread, the system was assumed to converge to the 
externally driven limit a = 1. That mechanism would, however, lead to a crossover 
to a = 1 with increasing At (Menezes and Barabasi 2004b). What is observed, is 
in fact not a crossover. There is no breakdown of scaling as a function of (/} for 
intermediate Afs as one would expect between the regime of the two limiting expo- 
nents (Menezes and Barabasi 2004b). On the other hand, it is well known (see, e.g., 
Gopikrishnan et al. 2000), that the number of trades Ni(t) is correlated. Individual 
trades tend to cluster together and this causes enhanced fluctuations in Ni(t). This 
mechanism sets in at time windows for which the probability for two trades to coin- 
cide is no longer negligible. The scaling law Q itself is preserved, but the exponent 
a is strongly affected. 

4 Conclusions 

In the above we have outlined a recent type of scaling analysis for the fluctuations of 
activity in complex systems. We have shown that systems can be classified according 
to the scaling exponents a. Then we have discussed how impact inhomogeneity and 
long range correlations give rise to non-trivial scaling exponents. Further research 
should clarify the interplay between fluctuations in the number of trades and in traded 
volumes/values in order to deepen the understanding of the market mechanism. 
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